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Abstract 

The exact solutions for the energy spectrum of the XX model with a periodic coupling and 
an external transverse magnetic field h are obtained. The diagonalization procedure is discussed, 
and analytical and numerical solutions are given. Using the solutions for period-two coupling, 
the free energy, entropy, and specific heat are calculated as functions of temperature and applied 
transverse external magnetic field. Their expressions show that below a particular value v and 
above a value u of the magnetic field \h\, the entropy and the specific heat vanish exponentially in 
the low temperature limit. 
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I. INTRODUCTION 



The physics of one dimensional chains provide a good example for the rigourous study of 
some properties of quantum magnetic systems. They are known to be useful for analysing 
many-body problems. Recently, real, model quasi-lD systems of atoms have become 
the object of experiments. From a experimental point of view, there are many quasi 
one dimensional compounds (such as the organic series (BCPTTF) 2 X with (X = AsF 6 , 
PF 6 ), the cuprate CuGeOs, or the (VO)2P2C>7 compounds, as well as other series as the 
TTFMS 4 C4(CF 3 )4 with M = Cu, Au, Pt, or Ni [1-7]). These are properly described as 
Peierls spin-chains. 

The experimental discovery of quasi-crystals [8, 9] has stimulated intense activity in the 
study of the periodic and aperiodic chains, with the aim of understanding their physical 
properties [10]. 

Recently, new, unusual properties in low-dimensional magnetic materials have been 
found [11-13] that can be explained in terms of many-body behaviour. It seems to be 
responsible for the occurrence of magnetization plateaus as a function of external magnetic 
fields [14, 15]. 

The procedure for the study of this type of structures was initiated by Leib et al. [16], 
by using the well known Jordan- Wigner transformation [17], and assuming that the chain 
is finite with set boundary conditions. 

One of the most important systems is the ID XY model (S = 1/2), introduced by Leib et 
al. [16]. It plays an important role in the description of many-body problems since it can 
be solved exactly. 

Two early examples of studies of ground state properties (e.g. susceptibility) for the XY 
model with alternating interaction are those by Perk et al. [18], and Taylor and Miiller [19]. 

Quantum critical phenomena in random XY chains have been studied using renormal- 
ization group methods [20], and by numerical methods [21, 22]. 

Recently, Derzhko et al. [23-25] have described the thermodynamical behaviour and prop- 
erties of the ground state of periodic non-uniform XY chains in a transverse field. Also, 
thermodynamic properties in one dimensional superlattices have been treated by de Lima 
et al. [26-28]. An excellent review is provided in Ref. [29]. 

As mentioned previously, systems of this type are interesting due to being strongly af- 
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fected by quantum fluctuations, and by their apparent simplicity. For example, low dimen- 
sional electronic materials are known to be very sensitive to structural distortions that are 
driven by electron-phonon interactions. These break the symmetry of the original ground 
state, which results in a new lower energy state where the electrons and ions are shifted from 
their original positions in a regular manner. This creates a periodic variation of the charge 
density that is called a charge density wave. This is the well know Peierls instability which 
opens a gap at the Fermi surface of the ID electronic chain, transforming a metal into a 
semiconductor. A similar effect is expected for ID spin chains, which are unstable against 
lattice vibrations. 

The problem considered here is a homogeneous chain where the structural distortions of 
the lattice are taken into account through a periodic coupling between a given site and its 
neighbouring sites. An external transverse magnetic field h is also applied. This permits us 
to study the diverse phases of the system. Thus, the Hamiltonian is, 

N N 

h = -J2^ + i(S!Sf +1 + sfsf +1 ) - hJ2s!, (l) 

i=l i=l 

where 

f2nm \ 

J M+ i = J i+0|i+0+1 = u + v cos I —fi-i J , (2) 

and the quantum spin operators are Sf = (1/2) (where a\ are the Pauli matrices at the 
site i), and a = N/m is an integer. 

This work is divided in four sections. In section II we begin with the Hamiltonian (1). 
The Jordan- Wigner transformation is applied, then the Hamiltonian is diagonalized with 
a superposition between creation operators. The same superposition is used between the 
corresponding destruction operators. This superposition is possible due to the Hamiltonian 
being isotropic, thus preserving the total transverse spin S z = Y^iSf. 

In section III, analytical results for the energy bands and gaps are given for the integer 
values a = 2 and 3. Graphs are provided that show numerical results for higher values of a. 
The degeneration of the bands into single levels is discussed. 

In section IV, thermodynamic quantities — free energy, entropy and specific heat — and 
their behaviour are discussed. The free energy per site is calculated as a function of tem- 
perature. The other thermodynamic properties are subsequently derived from this. Graphs 
showing the free energy and entropy are given. 
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An interesting result is found by taking the low temperature limit. It is shown that the 
entropy and the specific heat possess local maxima at \h\ = v and \h\ = u. They vanish 
exponentially with temperature for \h\ < v and \h\ > u. Both functions at the maxima vary 
as a square root of temperature, while in between they are proportional to temperature. 



II. THE MODEL 



The simplest solvable model (see for example Ref. [30]) is the isotropic spin 1/2 XY 
model. This model is based on a one dimensional lattice with a uniform interaction strength 
J between nearest neighbour sites, and uses periodic boundary conditions. As explained 
in the introduction, a correction to this Hamiltonian is provided by introducing a periodic 
strength, which simulates a vibration at the sites. Hence, the Hamiltonian from (1, 2) is 



N r 



H 



-E 



u + v cos 



(SfS 



x 

i+1 



qy qy \ h q z 

J i J i+1 j ) 



(3) 



where u is a uniform coupling and v is the amplitude of the periodic coupling. 

The spin operators Si act on a two dimensional space r]i = C 2 where we take the eigen- 
states of S 3 as a base. The total third component of the spin, 



N 

i=i 



(4) 



is a conserved quantity that can be used to describe the states of the system. 

The first task is to transform the Hamiltonian (3) into a fermionic Hamiltonian by means 
of the Jordan- Wigner transformation [17]. This is achieved by first defining the fermionic 
operators: 



at = K(l) , 
a\ = K(l) St, 



(5) 
(6) 



where Sf = (Sf ±iSf), and 



i-i 



i-i 



K(i) = e XP [i«T,s+ST\ = n(-2s;). 

V 3=1 J 3=1 

The new Hamiltonian and the conserved spin component S z become, 



A? 



# = -E 



n=l 



/27T 

u + v cos I — n 



^ (4a n+ i + 4 +1 a„) - h [ala n - J 



(7) 



(8) 



N / x 1 \ N 
S Z = Y, ( 4«n - = N f--2 J ( 9 ) 

n=l 

/ is the unity operator, and iV/ = S^Li (on^n) is the total fermionic number, which is also 
conserved. 

These a and at operators are true fermionic operators in the sense that |a|,a m | = 8i >m 
and |aj, aj^j- = = {ai,a m }. Furthermore, since 

of = 2a\ ai - 1, (10) 

the periodic boundary conditions require that 

cln+i = — Oi, q jv+i = — a i f° r states on which Nf is even; 

a N+i = Oi, a !v+i = a i f° r states on which Nf is odd. 
The next step is to transform the system to momentum space by performing a Fourier 
transformation using new fermionic operators bj and bj. These are related to the aj and a] 
by 

1 N 



h 3 = -fTf 12 eX P i ik U) l ] a h (11) 
1 N 

a i = -Trf 12 ex P l~ ik U) l \ b i, (I 2 ) 



and hermitian conjugate, where 



. !(,_£) with -o^M-dd, 

^ v 2/ ^ e = 1 if Afy is even, j = l,---,N. 

The new operators bj and 6] destroy or create a fermion with momentum k(j). 
In terms of these operators, the Hamiltonian (8) is 

N 

H = - U J2cos[k(j)} b]bj 

3=1 

- I E cos [k (j + f )] (e^ 6}^ + e — /" b) +m b 3 ) -h(^-12 fy^U) 

A particular case is the alternating chain that occurs when m = N/2. In this case, the 
Hamiltonian is 

N 

H = - u £>s[A;0')]&& 

TV /yy TV \ 

- i„ £ [km b]b j+N/2 -hl — -J2 h %\ ■ (15) 



The commutation rule of the Hamiltonian with the &j operator yields 
[H,b]] = {h -u cos [k(j)]}b] 



v 

- cos 
2 



e-^^6] +m , (16) 



Therefore, 6j gives the hermitian conjugate of this. As we can see, these commutation rules 
mix the operators bj and bj +m . Hence, we can redefine the index of bj in a new form, 

bj = bi t s with j = I + sm and I — 1, • • • , m and s — 0, • • • , a — 1. 

The periodic conditions are bi >a = fc^o, and fr^-i = fy, a -i- 
The sums in the index j will be transformed 



N 

E 



m a— 1 



and the functions k(j) 



EE witha =^ 

Z=l s=0 



2tts 



fcO') = k\l, s) = k(l) + 

(X 

Hence, the commutation rule (16) takes the form 
[H,b[ s ] = {/i-ucos^ 1 (/,*)]}&£, 



(17) 



(18) 



— - cos 
2 



fc 1 ^*)-* e^bl^ - V -cos k%s) + - e^b]^, (19) 



that we will write as 



with 



a-l 



r=0 



°s, s 


= /l 


— u cos fc 1 (Z, s) 






J 

c s-l,s 




V 

- COS 

2 


fc 1 (Z, s) - 


7T~ 

a. 


e iw/a 


J 

C s+l,s 




V 

- COS 

2 


k l (l, s) + 


7T" 

a. 




J 

C a-1,0 




f 

- cos 
2 


k\l, 0) - 


a. 


e «7r/a 


J 

C 0,a-1 




f 

- cos 
2 


k\l, 0) - 


a. 


e -i7r/a 


c * 

u r, s 


= 


otherwise. 







(20) 

For each /, these equations represent the elements of a hermitian a x a matrix that we 
write as 

C l = RJ- (21) 
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The eigenvalues of the these matrices yield the energy spectrum of the system. 

The case m = N/2 does not follow from (20) and must be calculated directly from the 
commutation relation with a = 2. The results for the elements of C l are 



-0,0 



h — u cos 



h + u cos 



k(l) 
~k(l) 



: 0j 1 = —iv sm 



k(l) 



The Hamiltonian is diagonalized using the operators 



a-l 



B 



i,p 



H S P, A 



s=0 



These are defined by the solutions of the eigenvalue equation 

[H, Bl}=EfBt p . 

The eigenvalues E l p are the eigenvalues of the matrix C l . 
Hence, the Hamiltonian becomes 



a— 1 to 



p=0 1=1 



(22) 



(23) 



(24) 



(25) 



Thus, the spectrum of this Hamiltonian is grouped into a bands, each of which has 
m = N/a levels. In the thermodynamic limit, keeping a fixed and finite, the sums over the 
index / become integrals. 



III. SOLUTIONS FOR a = 2 AND 3 

Case a — 2: 

For a = 2, the eigenvalues of C l are 

E? = h- yfv 2 + (u 2 - v 2 ) cos[A;(/)] 2 , 

E] = h + ^v 2 + (u 2 - v 2 ) cos[A;(/)] 2 . (26) 

The parameter / takes the values 1 < I < N/2; hence, when N — > oo, < k[l] < ir. 
Every value of the energy is doubly degenerate: E\ = E t N / 2 _ l . In Fig. 1 the energy of the 
two bands for the values h — 0, u — 1 and v = 0.5 is shown. For h — 0, the lower band has 
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Figure 1: Energy of the two levels Eq and i?i for u = 1 and v = 0.5 

i? G [— w, — f], while the upper band has E 1 G [f, u]. The difference in energy between the 
two bands is 2 min{w, v}. When u = v both energies are independent of / and every band 
collapses into a single level with values E n = =pw. 
Fig. 2 shows the difference, 

AE = El - Ei = 2^u 2 cos 2 k{l) + v 2 sin 2 k{l) 

The the minimum of AE represents difference in energy between the lower and upper bands. 




Figure 2: Energy difference between the two levels Eq and E\ 
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Case a = 3: 

In this case, the solutions for the eigenvalues of the C l are 



where 



El 
El 



48 V /( 2 ) J / ' 

^((-l+^v / 3)^ + (l+,v / 3)^), 

v (2) 



96 



V3)^ + (l-<V3)/«). 



./ 



(2) 



(27648 - 20736 + 6912 v 3 ) cos[3Jfe(Z)], 
-576 - 288 v 2 , 
l( / (D + ^ / (i) 2 + 4,(2)^ 1 



(27) 



(28) 



These coefficients are real. 

In Fig. 3, the energy levels of the three bands E°, E 1 , and E 2 (27) are plotted as a 
function of k(l) for four different values of v. The positions are marked on the x-axis where 
k(l) is, for I in the interval [l,iV/3]. 

The width of the bands are shown in Fig. (4) . They are given by the expressions 



-l-v + V9-6t;+ _M iiv>0 



1 +v 
~^2~ 



if v < 



(29) 




J k(l) 



Figure 3: Energy of the three bands E°, E 1 and E 2 for the values u = 1, and v = 0.1, 0.5, 1.0, and 
1.5 plotted as a function of k(l). 
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Figure 4: Energy band widths W(E°) = W(E 2 ) and W{E V ) 

It can be seen from Fig. (4) that with increasing v, the width of each band increases when v 
is negative and decreases when v is positive. The width is zero for v — — 1 and 2, meaning 
that at these two values each band degenerates into a single level. 

More generally, for an arbitrary value of a, the bands degenerate into single levels when 
v and u satisfy the condition 

1 



v 
u 



COS 



, j = a. 



(30) 



When this happens, the chain splits into a number of noninteracting parts. For example, 
j = 3 yields the solution v — — 1, and the chain is split into noninteracting subsystems with 
three sites each. If j = 1 or 2, then v — 2 and produces a chain where only the sites in the 
positions 3k and 3k + 1 (/c=integer) are interacting. 

The energy difference for three bands, AE = min \ Ef — E},\ = min \ E}„ — Eft,,], is given 

by 



AE 



3 + 3v - V9 - Qv + 3v 2 



(31) 



AE over the interval v — —1 to 2 is shown in Fig. 5. At fixed tt = 1 and large v, AE is 
proportional to v: 

AE u=ljV>>2 = Gv, (32) 
where the constant of proportionality 

3 - 



G 



(33) 
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Figure 5: Energy difference for three bands AE = min IE 1 — E°\ = min \E 2 — E l \ 
Case a > 3: 

Larger values of a are straightforward to deal with. Numerical results for the energy 
bands with a — 12, u — 10 and v = 3 are shown in Fig. 6. 

As follows from (30), that the bands can degenerate into a single level only when \v/u\ > 1. 
If we consider only v = ±u, then both cases satisfy this condition when a is even. However, 
if a is odd, this is true only for v = —u. In all cases that satisfy the degeneracy condition, 
the chain is split into noninteracting subsystems containing a sites. 



Energy 
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Figure 6: Energy Bands for it = 10, v = 3 and a = 12 
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IV. THERMODYNAMIC FUNCTIONS 



The partition function is defined by 

Z = tre-? H = Ee W (est|exp(-E P ,^f4 P ^ P ) |est> (34) 

| est) 

with p = 1/T, and T being the temperature. 

The most general state of the chain can be written as a tensorial product of the states 
of the sites which, in momentum space, are labelled by {I, p}. The state of every site, 
determined by the ket \o)i jP , can be occupied or empty: 

| est) = Hi iP \o)i tP , with o — 1, (occupied) o — 0, (empty) (35) 

Then, the trace is, 

]T (est| exp E P ,^f^ p B ljP ) |est) 

| est) 

= I] [<0| exp (-(3E*Bl p B hp ) |0) + (1| exp (pE^B^B^ |1>] (36) 

i, P 

Due to the boundary conditions, the function k(l) in Ef, in the product over l,p, is 
different when the number of occupied sites Nf in |est) is odd or even. However, when 
N — > oo the difference becomes negligible [16, 30], and 

Z = eW 2 n[l + exp(-/?£f)]. (37) 

This expression takes e = in fc(j) from (13) use it for Ef. 

The thermodynamic free energy per site / can be derived from the partition function as 
follows. 

/ = ~ ^ Z = -\ - E In [1 + exp (-f3Ef )] . (38) 

The thermodynamic limit is obtained by making iV — > oo and keeping a fixed and finite, 
whereupon the sums become integrals: 

^-V -> —V / dife. (39) 

When a = 2 the free energy can be calculated analytically. To perform the integrals for 
this case, we first change the site parameter / and define 

Ei = Ef = h- A(Z) for - j + 1 < I < — , 

N N N 
Ei = E) = h + X(l) for - - + KK — - and — + 1< I < — , (40) 

2 — — 4 4 — — 2 
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where 

A(Z) = \ju 2 cos 2 {k(l)) + v 2 sin 2 (A;(/)) = + (w 2 - v 2 ) cos 2 (A;(Z)). (41) 

Then, taking into account that E\ is an even function of/, the free energy in the continuum 
limit is 



/ = -\-—A r /2 \n(l + e-^ h - x ^) dk+ [ W In (l + e -W+*W>) rfA;" 



(42) 



Integrating by parts yields 

2 /3 V 2 " 1 + e-f('*+»> - V 1 " I J w 



where 



and 



h = L arCC ° S "V l + e-W+A) A (44) 



/■-* / /A 2 -v 2 \ e -^ +A ) „ 

/2 = y_ arccos (v^w j i + e -^ + A) rfA - (45) 

Both integrals can be performed by standard methods: see for example the appendix in 
Ref. [31]. Here, the function arccos is expanded as a series and integrated term by term. In 
the Ii integral, we must considerer the cases v > —h and v < —h. For J 2 , the corresponding 
cases are v > h and v < h. Then, Ii in the first case is 



h = T (| In (l + e"€) + £ T( n+ ^ M (n) i lfl (n)) , 



(46) 



where 



/2u 

01,1(0) 



Vm 2 — f 2 ' 
3w 2 + t; 2 



ai,i(2) 



6^/2v(u 2 -i; 2 ) 3 ' 
-5m 4 + 50uV + 3v 4 
80j2vHv? - w 2 ) 5 



7m 6 + 7m V + 301uV + 5v 6 . . 

3 = , 47 

448 ^2v 5 (w 2 - v 2 ) 7 



and 



u—v_ 



KM = J o YT^+t dx = " r 1 2 + n j PolyLog 1 2 + n ' J ' (48) 

with £ = (h + v)/T. 
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The expression for I\ when v < —h is more complicated. Following Refs. [32, 33], this can 
be written as follows. 



h 



arccos —\ 



lh 2 -v 2 
u 2 — V 2 



IT 



v - u E 

2 



u 



'h 2 — v 2 



arcsm 



h\ u 2 — v 2 



1 



+ 



where 



yj V - h 2 )(h 2 - V 2 ) 

\h\ 



oi, 2 (l) = 
ai,2(3) = 



+ ^ T n + i [i_(_i)«] a lt2 (n)i h2 (n), 

n=0 



h 



(49) 



u 2 - h 2 )(h 2 -v 2 ) 
h[2h 6 - Wh 2 u 2 v 2 + fr 4 (^ 2 + v 2 ) + 3mV(m 2 + t; 2 )] 
6^ - /i 2 ) 5 (/i 2 - t; 2 ) 5 



and 



00 



/o 1 + e x 
The first three elements of (51) are 

J2 



dx 



2 n _ 1 



r(n + i) C(n + i)- 



7T" . 77T 

11,2(1) = ?1 ' 2 ( 3 ) = 



4 

120' 



*1,2(5) 



317T 6 
^52" 



(50) 



(51) 



(52) 



Once again, \u — h\/T and \h — v\/T are assumed to be large. 

The I 2 integral when h < v can also be performed using the tables provided in Refs. [32, 
33], and a similar expansion to before. The result is 



I 2 = ( u -v)+ue(Ji-^ _r|ln(l + e^ 
-TET" + lfl 2)1 (n)i 2)1 (n), 



)/T 



(53) 



n=0 



where a 2j i(n) = — ai,i(n) from (47), and «2,i(^) = *i,i(^) from (48), with £ = (v — h)/T. 
The result for integral I 2 when h > v is 



/i arccos 



' /i 2 — f 2 
■u 2 — f 2 



+ W E 



arccos 



' h 2 — v 2 
■u 2 — 1> 2 



'1- 



1// 



+ E [1 - (-1)1 a 2 , 2 (n)i 2 , 2 (n), 



n=0 



where a 2 , 2 {n) = —ai t2 (n) from (50), and i 2>2 = i\ i2 from (51). 

The free energy is obtained by substituting I\ and I 2 into (43). Several examples are 
plotted in Fig. 7 for different values of T and v, and for u = 1, as functions of h. These 
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Figure 7: Free Energy plotted using u = 1 and v = 0.0, 0.25, 0.5, 0.75, 1. 



show that, as implied by the Hamiltonian (3), the free energy is an even function of h. This 
property is also shown explicitly by an alternative expression to (42) for the free energy: 

/ = -- / f ln(2cosh^ + 2cosh^) dk. 
n Jo \ T T J 

Notice that at low T, / becomes nearly constant between — v and v. 

The entropy per site is defined by 



(54) 



s = 



dT' 



(55) 



This can be obtained exactly from the expression for the free energy. To examine its be 
haviour we evaluate the partial derivative as follows. 

1 



7T 



L 



T 



where the function w(x) is 



w 



(x) = In (l + e~ x ) + 
15 



xe 



1 + e- 



(56) 



(57) 



This is an even function. Its magnitude is significant only when x is small: (w(5) = 0.04, 
w(10) = 0.0004). For large values of \x\, the expression simplifies to 

w(x) = \x\e~ lxl . (58) 

The integral can be evaluated in a similar manner to the one for the free energy. However, 
in this case the properties of the function w allow us to examine the behaviour of s in the 
low temperature limit. When \h\ < v , the integrand in (56) can be approximated by (58). 
Taking the most significant term, yields the following solution. 

s(h,T) T „ « I e -(--N)/ r _> 0. (59) 

7T 

A similar approximation can be made for \h\ > u: 

s(h,T) T „ w - e -^~ u ^ T -> 0. (60) 

7T 

To evaluate (56) for other values of h, we change the integration variable to 

h-X(x) h + X(x) 
V = j, , and y = (61) 

respectively, in the two integrals. Hence, the entropy is 

s = - / w{y) - dy 

71 V J(h-u)/T ^[ u 2 _ _ Ty )2] _ Ty )2 _ v 2] 

/• (/l+M)/T , x T ( T y -h) ,\ 

+ / w(y) V ; dy . 62 

^(h+t,)/T _ (/, _ Ty) 2 } [(h - Ty) 2 - v 2 } J 

This integral has local maxima at \h\ = v , where it can be approximated by 



S (| A | =v) = - f^ = 0.6475 . . (63) 

7r v 2(m^ — v l ) Jo ^fy y 2(w^ — ir) 

A similar approximation can be made when \h\ = u, where there is also a local maximum. 



Tu 

s{\h\=u) = 0.6475 ^ 2{u2 _ v2) (64) 

When v < \h\ < u — taking into account that w(x) has significant values around x ~ 0, 
and is rapidly decreasing otherwise — s can be approximated by 

1 T\h\ /■«> ^ A ^ ttT|/i| 

5 — 



/ «, y ) rf y = — 1 1 . 65 



K ^(u 2 -h 2 ){{h 2 -V 2 ) J-oo ^ [u2 - h 2 ){h 2 - V 2 )' 
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Here, we can use 



-f$rfx = 0.6475, - r w(x) dx = - 

71 JO Jx 7T J-oo 3 



Thus, combining all of the above, the entropy is 

I e -(v-\h\)/T 

0.6475 



s(h,T) T ^ 



Tv 
71 ZT\ 



p(u 2 -v 2 ) 
nT\h\ 



3xJ(u 2 -h 2 )(h 2 - 



0.6475 



/ Tu 
2(u 2 - v 2 ) 



7Y 



<\h\-u)/T 



if \h\ < v 
if \h\ = v 

if u > \h\ > v 

if\h\ =u 
if \h\ > u 



(66) 



(67) 



The entropy is plotted using four different values of v, each at three different temperatures 
in Fig. 8. Fig. 9 shows logarithmic graphs of the temperature dependence of the entropy 
when u — 1 and v = 0.5, for five values of values h. 

The specific heat per site is 



C(h,T) T ^ = T 



ds 
df 



v — 



v T 
0.3238 



-{v-\h\)/T 



3< 

0.3238 



/ Tv 


V 2(« a " 


v 2 ) 


7TT\h\ 




-h 2 ){h 2 


-v 2 ) 


1 Tu 


V 2 (« 2 - 


v 2 ) 



-« p -(\h\-u)/T 



if \h\ < v 

if \h\ = v 

if \h\ > v 

if\h\ =v 

if \h\ > u 



(68) 



It can be seen that both functions (the entropy and the specific heat) decay exponentally 
as e~ 1//T to zero when T — » 0, for \h\ < v and \h\ > u. This also shows that the internal 
parameters of the system, v and u, can be derived by measuring the specific heat as a 
function of h. 

Conclusions for other values of periodicity a can be generalised from this result. The 
specific heat will be exponentially small for values of h between the several energy bands. 
Thus, measurements will provide information about the periodicity and its amplitude. 
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